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The crossover from thermal hopping to quantum tunneling is studied. We show that the de- 
cay rate F with dissipation can accurately be determined near the crossover temperature. Besides 
considering the Wentzel-Kramers-Brillouin (WKB) exponent, we also calculate contribution of the 
fluctuation modes around the saddle point and give an extended account of a previous study of 
crossover region. We deal with two dangerous fluctuation modes whose contribution can't be cal- 
culated by the steepest descent method and show that higher order couplings between the two 
■ dangerous modes need to be taken into considerations. At last the crossover from thermal hopping 

I , to quantum tunneling in the molecular magnet Mni2Ac is studied. 

<g ■ PACS numbcr(s): 75.45. +j, 75.50.Tt. 
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I. INTRODUCTION 



The decay of metastable states in macroscopic systems is a fundamental problem in many areas of physics, such 
as macroscopic quantum tunneling in Josephson system Jl]J^] , violation of baryon-lepton in Weinberg-Salam model, 
nucleation in first order phase transition theory J^^J] and more recent magnetic quantum resonant tunneling || . The 
crossover from thermal hoping to quantum tunneling has been studied intensively. Using functional integral approach, 
Affleck first demonstrated the transition can be found between classical regime and quantum regime |2g| ] . Larkin and 
Ovchimikov also suggested it and gave a formula determining the boundary of first order transition and second order 
one |?],|8|]. Grabert and Weiss discussed the phase transition in the presence of dissipative effects of the environment 
m some detail §-0. 

At high temperature, the decay of the metastable state is determined by process of thermal activation, which is 
governed by the Arrhenius law, ^exp(— AU/T) where Wo = \/ — U"(xq)/M is the well frequency and AU is the 
, barrier height. While at T = the particle can escape from the metastable state due to quantum tunneling, the rate 
of which goes as exp(— B) where B is the WKB exponent. Ignoring the prefactor and equating the exponents, one 
0"S \ obtains the estimate 

G\ . 

T (0) = ^, (1) 



AU 



where the superscript at T means that the ground state tunneling is considered. For T > lg , one has practically 
r = Ttherm(T), whereas below the transition T = T quan is independent of temperature. 

It turns out that for common metastable or double-well potentials, such as cubic or quartic parabola, below the 
crossover temperature To the particles cross the barrier at the most favorable energy level E(T) which goes down 
• • ■ from the top of the barrier to the bottom of U (x) with lowering temperature. The second-order transition from the 
classical thermal activation to thermally assisted tunneling (TAT) is smooth and the transition temperature is given 
by 



X 

& ■ rf = f£ (2) 



where u>b = \/—U"(xb)/M is the barrier frequency, and Xb corresponds to the top(the saddle point) of the barrier. 

In ref. ( |12j ) , Chudnovsky stressed the analogy of this kind of transition phenomena with ordinary phase transitions 
and analyzed the general conditions for both types of quantum-classical transitions. For the second-order transition 
the period of oscillation t(E) in the inverted potential —U{x) increases monotonically with energy E from the top 
of the barrier. If t{E) is not monotonic, the first-order transition occurs. The escape rate can be conveniently 
represented in terms of the effective temperature defined by 
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r ~ exp[--^-] = CXp^S'min/S,]. (3) 

The actual dependence of S m i n (T) goes along the minimum of these two actions ( sphelaron and periodic instantons 
) and the first-order transition occurs at T — The first derivative of S m i n (T) is discontinuous at providing 
that the crossover from the thermal to the quantum regime is the first order transition. 

The second-order transitions are common, whereas the first-order ones are exotic and have to be specially looked 
for. Nevertheless, a number of systems and processes that show first-order transitions are already known, e.g., a 
SQUID with two Josephson junctions Jl3| |, false vacuum decay in field theories P,|l4|-|l6|| and depinning of a massive 
string from a linear defect All these systems have more degrees of freedom than just a particle, thus search 

for a physical system equivalent to a particle in a potential U(x) leading to the first-order transition of the escape 
rate seems quite actual. Qualitatively it is clear, how U(x) look like: The potential should change slowly near the top 
and the bottom, but is rather steep in the middle. In this case, as for the rectangular barrier, tunneling just below 
the top of the barrier is unfavorable, the TAT mechanism is suppressed, and the thermal activation competes with 
the ground state tunneling directly, leading to the first order transition. 

Quantum tunneling of the magnetization(QTM) has become a focus of interest in physics and chemistry because 
it can provide a signature of quantum mechanical behavior in macroscopic system |3,[l2],^8 19-24|. At low enough 
temperature, it has been demonstrated that the vector of the magnetization formed by a large number of spins in 
magnetic system can coherently tunnel between the degenerate minima of magnetic energy. Theoretical suggestions 
have led to a number of experiments which seem to support the idea of magnetic tunneling. Since the M1112 Ac complex 
magnetic molecule provides a more suitable model for the magnetic quantum tunneling, extensive works have been 
performed to demonstrate the QTM in large spin molecules po] p7{ . On the other hand, Mn^Ac molecular is one of 
the very few examples which could exhibit the first-order transition. 

We derived a compact formula for decay rate which is valid for the entire range of parameters of the interest in the 
problem of MQT. The quantum-classical transitions of the escape rates in the dissipation systems are investigated by 
the periodic instanton method. Applying the periodic instanton method, we showed that the first-order transitions 
occur below the critical external magnetic field h x — j for Mn^Ac molecular which is in good agreement with earlier 
works 
and ca 



28|.|23| . The results of the application of a previous method is developed for dealing with the quasi-zero modes 
culating of decay rate in the crossover region which is beyond the steepest descent method || |ll| . In the 
crossover region of second order transition h x > j, we take fourth-order terms into account to include nonlinear 
couplings between the modes and obtain the universal law in this region. The point h x = -7 is the boundary of first- 
order and second order. At this point the fourth order couplings between modes disappear and sixth order terms must 
be considered. For h x < \ before the eigenvalues of the two quasi-zero modes reach zero, the first order transition 

occurs and there is no universal law in common. While T (1) is not far from T (2) ( such as T (1) = 1.078T (2) for h x = 0.1 
), the two dangerous modes also play important roles on the tunneling rate and need to be calculated carefully. 



II. DECAY RATE IN THE CROSSOVER REGION OF SECOND ORDER TRANSITION 

The partition function can be written as a functional integral over periodic paths where the path probability is 
weighted according to the Euclidean action 

S= / gIt[-Mx 2 + U(x)} (4) 
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where k (r) = j^m X) £ i v n) ex P {i v n r ) 1 v n — litn/fth, and £ (i/„) = 7 (i/ n ) \ v n | is related to the frequency- 



n— — 00 



dependent damping coefficient 7 (u n ) . U(x) is a metastable potential with a local minimum at x = and a local 
maximun at 2; = Xb- We use lor to denote the solution of the following equation ui R + uir^(u!r) = ui^ where 
cob = — U" (xb) /M characterizes the width of the parabolic top of the well. In classical limit, i — ► 00 , the steepest 
descents method is available: 

SS [x (r)] =0, x(0)=x (f3h) . (5) 
The fluctuation modes about the saddle point are expanded using ^! n , x — x c (r) + ^ Y n ^ n , where Y n are fluctuation 
amplitudes and ^„ are modes of the spectrum : 
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9 n + U [x c (r)] * 

* (/3^) = # (0) 



(6) 



j^ImF. Above the crossover 



According to the metastable-decay theory quantum tunneling rate has the form Y - 
temperature T c , the decay process comes from the thermal activation V = jr-£-ImF where (3 C — 2it/uir p8| . In 
ordinary case, the one loop correction which results in a prefactor of the WKB leading order exponential does not 
enhance the tunneling significantly and the transition rate is dominated by the WKB leading order exponential. Near 
transition point the imaginary part of the free energy has a common form : 



I m F 



1 



,(0)2 



A 



-/ c exp 



Sr. 



where S c is just the WKB leading order exponent, uo n = ujq + v\ + u n j (i/ n ) , u, 
comes from the two quasi-zero modes which need to be calcaulated carefully and 

f c [UQ,U b ] = ] | 



(6)2 



n=2 



(b) 



(7) 



-u% + vl + v nl (u n ) , 1/A 



(8) 
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= const 
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where A^ = — ^ ± 



\± — 7 4- 



and Wo = \J —V" (0) /M. Years ago Grabert and Weiss discussed 

the transition rate in the presence of dissipative effects of the environment in some detail 1-0. They found an 
unstable mode besides the zero mode near transition point and calculated it carefully |]6|-|l0|] . Near phase transition 
point the fluctuation modes about the saddle points include two dangerous modes whose contribution can't be calculated 
by the steepest descent method and it is necessary to consider higher order couplings between the two dangerous modes 



(2) 

A. Beyond steepest descent for T > T 

(2) 

Above T , the decay process is dominated by the saddle point called sphelaron x = x b . Considering the fluctuation 
modes around it, we have the periodic paths near the saddle point g-Q 

2tt 2ir 

x = x b + Y +Yl 1 -v / 2 sin— t + Y x V2 cos— t + ... (9) 

Lj Li 

r 2im r 2nn 

+y_„V2sin — — t + Y n \J2 cos — — r.... 

There is a mode with negative eigenvalue to^ 2 — —u>% = —U"/M which is the key mode giving contribution to the 
imaginary part of the free energy. From the steepest descent method the result of the partition function of sphelaron 
solution is written into 

Z= J D[x(r)]exp{^lH} (10) 

x(0)=x(/3h)=0 



= N [ l[dY n exp{^l] 

11 N -s c 

n#0 



(0)2 

1*0 



where N — fiulQ |p9|] . The eigenvalues of the two lowest positive modes are Ai = to^ 2 — uj^\ 2 = v\ 

Second order transition occurs when the eigenvalue of the lowest modes is equal to zero Ai = as temperature 



R- 
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decreases, so it is defined that: T c = <jJr/2it. Near the transition point, the eigenvalue of the lowest positive modes is 



Ai 



-as 



(11) 



where e = (1 — i^-) and a — uj\ + lo 2 r [1 + dj (uj r ) /du} B ] . 

To regularize the divergent integral we have to add terms of fourth order in the amplitudes Y±\. After expanding 
the potential about the barrier top 



E 



CiX 



(12) 



where Cj = (x — x^) / i\ , we obtains the action 
S[q] 



£ At 1 {± c (r) + Y n V^J + V L c (r) + £ Y n %, 



(13) 



h(3 AU + -h(3m 



AS 



where 



As = pn 



(14) 



B t = |c 4 + tj^t - TrPrm and ^2 



<)(:- 



(6)2 



4^ 2 -w| 



3wfj. After integrating the amplitudes Yq and Y n , we deal with the 

2 



2^ 1 2M^f 4M W < 6)2 

quasi-zero modes to considerate the fourth term : B4 (Yf + Y^Y ■ Introducing the polar coordination : pcosd = Y\, 
psinO — Y-i, we get 



1 = JL 

A 2tt 

KH/7T 



2*4 



erfc (— «e) exp (k 2 £ 2 ) 



(15) 



where k 



B4 



transition and the second order one in ref. ( |§). For S4 < the integration in eq. ( |l5| ) is divergent and this kind of 
divergence will be discussed in first order transition cases. It is obvious that the dissipation may change the boundary 
between the first order transition and the second order one. Then we have the transition rate from eq. (M) and eq. 

(i 



!)<•- 



has been defined as the boundary between the first order 



%(°' 2 ^ erfc (-ks) exp ( K 2 e 2 ) f c [c , u> b ] e~ u ^ T . 



(16) 



Away from the crossover region the result goes back to the classical decay rate u / ksT . 



B. Beyond steepest descent for T < T 

Below the crossover temperature, the saddle point is named by periodic instanton or thermon. There are also 

(2) 

two dangerous modes about this saddle point near T : one is quasi-zero mode which is associated with a phase 
fluctuations of the periodic instantons with the eigenvalue and eigenstates of <*4 = 2ae and ^2 = \/2 sin {oj^t); 
the other represents amplitudes fluctuation and gives large contribution to partition function with uj^ 2 = and 
cos (w&r). The quasi-zero mode just takes place of "soft mode" which restores symmetry and the zero mode 
of Goldstone mode which reflects the freedom of phase. This is just the character of Global U(l) symmetry broken. 
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(2) 

Near T , this kind of classical periodic trajectory of thermon may be written as a Fourier series 



(r) = ^ [X n cos (i/„t) + X^ n sin (y n r) 



(17) 



The periodic paths near the saddle point are similar to that of eq. (^) x — x c (r) + E ^n^n — E (Xn + X n ) \t n . We 

n n 

define the amplitudes into another form = Y n + X n and the action is 



S [q] = h[3 A U + -Tiftm 



n— — oo,n^±l 



AS 



where 



AS = -/3hmLu[ b}2 {Y{f + -ptimJ^l 2 (Y^f 



(18) 



(19) 



-Ba 



(Y>f + (Yuy 



In term of Y±-y we obtain the tunneling rate below T 



(2) 



1 (0)2 
F _ 1 ^0 , r j -AU/k B T 

Tip Wb A 



(20) 



where 



1_ _ k^/tx 
A ~~ 2^1" 



erfc (— ke) exp (k e ) and k 



, y -jf- ■ In term of k the size of crossover region is defined 



T w_ T 



< — in which the crossover occurs. It is obvious that there is symmetry 

(2) 



T—TX 



Ti '-T . 



M7T 



in the crossover 



region | -J-pj- |< - 



T (2) _ T 



Away from the crossover region | (2) | > i , the tunneling rate reduces to the standard form 



Z b = 



E -^ )2 
i V 

A / v (6^2 ' 
n 5^0,1 



(21) 



where-^ = y ^%(/3^-) which is known by Faddeev-Popov technique. We can reach it only through normalizing the 
eigenfunction of zero mode while now N may be different from the upper result for all modes being normalized 



N 



The concrete parameter 



E »i b] 



can be calculated only when T -> T (2) @M|| 



C. Universal law in crossover region 

Beyond steepest descent method, we have the formula eq. (|7|) and eq. (||), which is only needed in the crossover 
region 



T _ T (2) |< r W/K 



( 2 ), 



(22) 



where k = M ^ R yl ;g> 1. It has been pointed out that there is a universal law in the crossover region of second order 



transition pJlQ|| . We use the following quantity to show the universal law 

y = Tcxp(AU/k B T). 



(23) 
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which is a function of e but independent of the temperature T. According to the formula eq. (0) and eq. (|8|), we have 
the universal law in the scaling region 

v/Vo = F (£/£„) (24) 
where F (0 = erfc (0 exp (£ 2 ) , £ = T — T (2) , £o = T (2) /« and y = ^i^±f!lJK fc [w0j Wfc] . 



III. DECAY RATE IN THE CROSSOVER REGION AT THE BOUNDARY OF SECOND ORDER AND 

FIRST ORDER TRANSITION 



£>4 = is the boundary of the second order transition and the first order one and the potential looks different from 
that of £?4 ^ : The potentials changes slowly near the top and the bottom, but are rather steep in the middle. 
Because there is divergence of n — ► oo at the point of B^ — 0, the formula ( Jlq ) isn't available. 

(2) 

Above T , the decay process arises from sphelaron x = Xb, too. So we have the same periodic path near the saddle 
point as formula (||). but the interactions terms of modes differ from eq. (|l4|). We consider the sixth order terms of 
the two dangerous modes 



AS = (3H 



Muj\ 



-B 4 {Y? + Y 2 1 ) 2 + B 6 (Y 2 



Y 2 i)3 



(25) 



where B e = |c 6 



Mu 



Near the two-phase point not only the second term £ d£ )2 Y% tends to 



n=±l 



(Remember here — > 0, w/ien temperature turns to T c ) but also B4 is a small quantity. We consider the sixth 
order terms of the two dangerous modes Bg(Y 2 + Y 2 1 ) 3 . At the point B4 = 0, there is no fourth term and the formula 
is reduced to 



1 

A 



[/3B 6 }' 1/3 / (ft exp [- (t 3 - 3k' et)] . 
Jo 



where k' = [(3B 6 ]~ 1/3 . Below T c , we transform Y n to Y n + X n and have the same form of 4- 

The universal law in the crossover region | T — T c |< T c /k! is also defined as y/ya 
/ °° dt exp [- (t 3 + 3#)] , C = 7 1 - T c , - T c /k' and 



■1/0 



w / 2 
12^6 ^ " 

[/3S 6 ]" 1/3 : 



(26) 

Fii,';,,) where F (0 = 
(27) 



, r (a;) is Gamma function. 



IV. THE EFFECT OF MODES OF THE FIRST ORDER TRANSITION B 4 < 



The actual transition occur at the temperature when the two saddle points have the same action. For B4 > the 

particles tunnel through the barrier at the most favorable energy level E(T) which goes down continuously from the 

top of the barrier to the bottom of the potential with lowering temperature. This corresponds to the second-order 

transition from thermal activation to thermally assisted tunneling ( TAT) with no discontinuity of dT/dt at To. And 

(2) 

the transition temperature is given by T . For the cases of B4 < 0, just below the top of the barrier is unfavorable 
for tunneling, the TAT mechanism is partially suppressed and the first order transition occurs. Because the top of the 
barrier is wider for B4 < 0, the particle doesn't have more advantage to tunnel through the barrier from the higher 
excited states than the lower ones. 

Let us discuss the dangerous modes near the first order transition point: Before the eigenvalues of the two quasi-zero 
modes reach zero, the crossover occurs and tunneling process is dominated by periodic instantons just below the top 
of the barrier. In this case, there is really no universal law. Above T c which is near lor/2tt, the sphelaron's two 
dangerous modes may also play important role on the tunneling rate at the transition point. Above the crossover 
temperature, the quasi-zero modes of sphelaron have the same eigenfunction as other cases (B 4 > ) : V2 sin - j-r 



G 



andv^cos =£t. While there are many distinguished differences of the interaction of modes between the case of B4 > 
and that of B4 < : The interaction between modes is attractive of the lower order terms. 

After integrating the amplitudes Yq and Y n , we deal with the quasi-zero modes to considerate the forth term 



1 = JL 

A 2tt 



x exp —/3 



dYidF_i 
1 



(28) 



2 Tnwl b)J Y? + ^mo/T^-i + B A (if + Y* )' 



a/7T 2 er ^ (~ K£ ) ex P ( k2£2 ) 
2cJ 



where k 



A/u 



-J^- and e = (1 ^57) . The Quasi-zero modes give a divergent contribution to the partition function 

00 00 c 

which can be distorted into complex-plane : J exp (x 2 ) dx — J exp (a; 2 ) dx — J exp (a; 2 ) da; = i}£ZL — erfi (c). Here 

c 

c 

erfi (c) = J exp (a; 2 ) da;. Because the real part of the function erfi (— ne) goes to zero as e — > 0, a dramatic phenomenon 


arises - first order transition may suppress quantum tunneling rate (without considering higher order terms )! 

Because the periodic trajectories of periodic instanton can't be written as a Fourier series if they are far below the 

.(i) 



top of the barrier, the formula (|28|) is useful only above the crossover temperature T 

Now we consider a weak first order phase transition which is under the condition T /Tq 2 -* < 1 — For a weak 

first order phase transition, the crossover temperature Tq 1 ^ is in the region (T — Tq 2 ^) < Tq / K which is shown in 
Fig. 3. We may define a remaining crossover region above the crossover temperature 



T (1) <T< T (2) +T (2) /k- 
We also use the quantity y = T exp (AU /ksT) to show the universal law as 

y/yo = 

("S+"g) fW* r,... ,..1 T7_ <rW /t( 2 ) 



where F(£) = erfi(£) exp(£ 2 ) and y = ^ 2 N,^h]- For the cases T^> /T^> > 1 - 1/k, the first order 

transition is rather sharp, there is no crossover region at all. 



(29) 



(30) 



V. SCALING LAW OF QUANTUM-CLASSICAL TRANSITION OF THE ESCAPE RATE IN MN12AC 



A rather simple and experimentally important system which may exhibit the first-order transition is the uniaxial 
spin model in a field parallel to x-axis H x described by the Hamiltonian 



H = -DSi - H X S X 



(31) 



which is generic for problems of spin tunneling. This model is believed to be a good approximation for the molecular 
magnet Mn^Ac with D the anisotropy constant. Using the method of mapping the spin problem onto a particle 

_ p 2 



(32) 



problem [|22|j25|], we have the equivalent particle Hamiltonian as H = ^ + U{x), where 

U(x) = S 2 D(h x coshx- l) 2 



and m = ^j, h x — and S> 1. The minimum of the effective potential, xo = cosh 1 has been moved to zero. 
Integrating the equation of motion with imaginary time variable one obtains 



x p( t ) — 2 tanh 1 



tanhxi sn ( -, k 



(33) 



tanh x\ 2 



1 - K T VE 7 



(34) 
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where sn (r, k) is the Jacobi elliptic function with modulus k and the complementary modulus k' = \fl — k 2 , E' = ^f^y . 
£p is the characteristic length of the periodic instanton determined by the following equation 

l-hl 

^ = 52^2! v f ^+k 2 )hi-k' 2 (i-hi)+h^ihik*+k<* y ] . (35) 



(l + k 2 )hl + h x ^/4hlk 2 + k' 4 



This description corresponds to the movement of a pseudo-particle in the inverted potential —U(x) with energy — E. 
The periodicity of the solution (^) is t(E) = 4K/3, K/3 — K(fc)£p, where K(fc) is the complete elliptic integral of the 
first kind. The Euclidean action of the periodic instanton configuration in the whole period is 



Sp 



- 2 
m x p 

-P 



dr = W (36) 



where 



W = —^[(a 4 -k 2 )Tl(a 2 ,k) + 

fc 2 K(fc)+a 2 (K(fc)-E(fc))] (37) 

where a 2 = tanh 2 x\ < k 2 and E(fc), n(a 2 , k) are the complete elliptic integral of second and third kind, respectively. 

The period of oscillation t(E) in the inverted potential — U{x), i.e., the periodicity of the periodic instanton solution 
(|33|), t(E) = 4K(fc)£p, can be equivalently calculated by 

dW(E) r— f Xl dx 

T ^ = JW 1 = ^ / 7777 x 5 " 38 

dE J_ xt y/U(x) - E 

Near the top of the barrier, k — > 0,K — > 7r/2, the equation ( |3^ ) yields the previously known result r = 2-k/lui,. Near 
the bottom, one has k — > 1, and r logarithmically diverges. For h x < 1/4, the dependence r(E) is non-monotonic, 
and the transition is first order. 

Since both Sq and S p are assumed to be large compared to h, the smaller of the two determines the actual escape 
rate. The calculation of the temperature dependence of iSmin is depicted in Fig.l for h x = 0.1. The solid line 
corresponds to the periodic instanton action S p while the dashed one corresponds to the thermodynamic action Sq. 
The actual dependence of S'min(T') goes along the minimum of these two actions and the first-order transition occurs 
at T = T (1) satisfying T 2 < T (1) < T u where k B T x = %/t{E 1 ), k B T 2 = H/t(AU). The first derivative of S min (T) is 
discontinuous at Tq \ , providing that the crossover from the thermal to the quantum regime is the first order transition 
on temperature. Quite recently an effective free energy F = acf> 2 + bej) 4 + c</> 6 for the transitions of a spin system was 
introduced |24| as in the Landau model of second order phase transition. Here a changes sign at the phase transition 
temperature T = Tq 2 ' = y/ h x (l — h x ) and 6 = corresponds to the boundary between first- and second-order 
transitions. There indeed exists a phase boundary between the first- and second-order transitions, i.e., h x = j, at 

which the factor in front of <fi 2 changes the sign. At h x = 0.3 the minimum of F remains AU for all T > T (2) . Below 
(2) 

T it continuously shifts from the top to the bottom of the potential as temperature is lowered. This corresponds to 
the second order transition from thermal activation to TAT. At h x = 0.1, however, there can be one or two minima of 
F, depending on the temperature. The crossover between classical and quantum regimes occurs when the two minima 
have the same free energy, which for h x — 0.1 takes place at Tq — 1.0787g . 
Near the top of the barrier the potential has the form as 

U(x) = S 2 D(h x coshx- l) 2 (39) 
= 2S 2 Dh x - S 2 Dh x (h x - l)x 2 

+C42; 4 + CqX 6 ... 

where C4 = S Dh * 1 ^ 1 * — il and c@ = 2S 2 Dh x (h x — ^)/45. The potentials changes slowly near the top and the bottom, 
but are rather steep in the middle for the cases C4 < 0. 

h x = 1/4 is the boundary of the second order transition and the first order one. When h x > 1/4, we have the 
formula eq. (|2^) and eq. (pi|), which is only needed in the crossover region 

| T _ Tq (2) |<t (2) /k (40) 
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where 



K = [h x (l - h x 



,3/4 



2ttS 



h x (h x - i) 



S 1 / 2 » 1. 



(41) 



It has been pointed out that there is a universal law in the crossover region of second order transition [tUj , !!] . According 
to the formula eq. (P3|) and eq. (M), we have the universal law in the scaling region 



y/yo = f (e/^o) 

where F (£) - erfc (£) exp (£ 2 ) , £ = T - T (2) , £ = T (2) /« and y 



(42) 



2tt 2 
(2) |^ ™(2 



For the case h x = 1/4, we have the crossover region | T — Tq |< Tq j /V where 



k = 6 



7171a; (1 - 7l x )5 



-i 2/3 r 



h x(K - — ) 



-1/3 



5 2/3 > 1. 



(43) 



^From the formula (p7|), we have the universal law y/yo = -^(C/Co) in the crossover region where F (£) 
/ °° cftexp [- (r! 3 + 3ft)] , £ = T- T (2) , £ = T (2) /«' and 



24tt 



-1/3 



r ( g ) fc [wo,w b ] , 



r (x) is Gamma function which is shown in fig. 2. 



Now we consider a weak first order phase transition which is under the condition Tq 1 ^ /T < 1 



(44) 



For 



a weak first order phase transition, the crossover temperature Tq is in the region (T-T (2) ) < T (2) /k which is shown 



[Mi-Ml 3/4 

f)<T «/ 

in Fig. 3. We may define a remaining crossover region above the crossover temperature 

T (1) < T < T (2) + T (2) /k. (45) 
We also use the quantity y = rexp (AJJ/fe^T) to show the universal law as y/yo = -FXC/^o) where -F(f) = 

,(1) /rr(2) 

is rather sharp, there is no crossover region at all 



/ 2 2 n / h»e>»-j) 

erfi(f)exp(<f) andy = l^^^/ch,^. For the cases T (1) /T (2) > 1- V the first order transition 



VI. CONCLUSION 

In this paper we have shown that the decay rate T can accurately be determined near the crossover temperature 
in dissipative systems. Besides considering the WKB exponential, we have calculated contribution of the fluctuation 
modes around the saddle point and have given an extended account of the previous study of crossover region [p|-pd| . 
Near phase transition point the fluctuation modes about the saddle points include two dangerous modes whose 
contribution can't be calculated by the steepest descent method and the higher order couplings arc considered between 
the two dangerous modes, near the point of 



B, = -c 4 



_9c|_ 
2Mw$ 



9cj 



AMuj. 



(6)2 



= 



(46) 



sixth order need to be considered. The results can be easily used in Mn^Ac of which the equation B4 = 
S(S+i),D6 x (fe x — ^2 _ q gj veg phase boundary point h x = j in good agreement with earlier works. 

Another example is biaxial anisotropic ferromagnetic model H = K\S 2 + K%Sy which describes XOY easy plane 
anisotropy and an easy axis along the x direction with the anisotropy constants K\ > K2 > 0. Mapped onto a particle 
problem, the equivalent particle Hamiltonian is 



H = 



AKi 



x -K 2 S(S+ l)sn 2 (x,X) 



(47) 







where sn (r, A) is the Jacobi elliptic function with modulus A = K2/K1 . From the equation B4 = g2g ( s+1 )(i x )( x 2A ) = 
0, we obtains the phase boundary point A = | whichs confirm the results in ref. ( |02| ), 

The results of the application of a previous method 11 1 is developed for dealing with the quasi-zero modes and 



calculation of decay rate in the crossover region which is beyond the steepest descent method. The decay rate is valid 
for the entire interesting range of parameters in the problem of MQT. 
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Figure Captions 

Fig. 1 

First-order transition from the thermal to the quantum region for Mni2 molecular: h x = 0.1. 
Fig. 2 

Scale factor F (£) at different sides of phase transition point of h x = 1/4. 
Fig. 3 

The remaining scaling region of a weak first order phase transition is T <T< + / k. 
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